Is a step in the primordial spectral index favoured by CMB data ? 
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A sudden small change in the second derivative of the inflaton potential can result in a universal 
local feature in the spectrum ofprimordial perturbations generated during inflation. The exact 
solution describing this feature [y is characterized by a step in the spectral index modulated by 
characteristic oscillations and results in a large running of the spectral index localized over a few 
e-folds of scale. In this paper we confront this step-like feature with the 5 year WMAP results and 
demonstrate that it provides a better fit to this data than a featureless initial spectrum. If such a 
feature exists at all, then it should lie at sufficiently large scales fco < 0.003Mpc~^ corresponding 
to I < 40. The sign of the effect is shown to correspond to the negative running of Us localized 
near this scale. This feature could arise as a result of a 'mini-waterfall'-type fast second order phase 
^~>' transition experienced by an auxiliary heavy field during infiation, in a model similar to hybrid 

infiation (though for a different choice of parameters). If this is the case, then the auxiliary field 
should be positively coupled to the infiaton. 
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INTRODUCTION 



i-^ ' The great precision of current cosmological data and the enormous volume of data expected in coming years leads 

"T^l to the hope that cosmological parameters will soon be determined to great accuracy. Since cosmological parameters 
Q ■ are intimately linked to an underlying theoretical model, the increasing depth and sophistication of cosmological data 
sets, especially those associated with the Cosmic Microwave Background (CMB), will, in all likelihood, lead to a more 
^ t refined and deeper understanding of such important issues as the form of the initial perturbation spectrum, the energy 
^1 ' scale of inflation, etc. 

Indeed, the cosmological world view which has emerged during the past two decades has been quite intriguing! New 
04 ' cosmological data has, while providing support to earlier constructs such as inflation, also paved the way for radically 
^ , new propositions such as dark energy. 

^T ■ An important issue concerns the form of the primordial spectral index ns{k) = dlnP{k)/d\nk. The large quantity 

i_: I of CMB data which has become available over the past few years indicates that the departure of the spectral index 
fK^ . from exact scale invariance is likely to be small, \ns{k) — 1| <C 1, which is in good agreement with predictions of the 
simplest inflationary scenarios (see e.g. [3| for a recent model-independent analysis of observational data assuming 
local analytic behaviour of the inflaton potential V{(f>) in the observable window). An associated question concerns 
the value of the running of the spectral index a{k) = dns{k)/dlnk. Most inflationary models with smooth inflaton 
(^ ] potentials do predict a non-vanishing value of a{k), usually |a(fc)| ~ \ns{k) — Ip ^ 1. However, recent WMAP^ 
; • ■ results [3, 3, [1^ may be suggesting a somewhat larger value |a(A:) | ~ |ns (fc) — 1 1 . Taken together with other 'anomalies' 
such as the 'Archeops feature' at / ~ 40 [3, 0, Q, these recent data may be providing a subtle hint that inflationary 
models are slightly more complex than the simplest single-field models suggested during the early 1980's. 
$H Note that the only way for a large value of the running a(k) to be accommodated within the inflationary paradigm is 

for it to be localized (i.e. have the form of a feature), since, otherwise, the number of inflationary e-foldings turn out to 
be too small [6|. Generation of features in the primordial inflationary spectrum has been the subject of considerable 
discussion. A partial list of references may be found in [3, M, M, [l3, IHI- Several model independent attempts to 
search for such features in the primordial spectrum directly from the WMAP data have led to some positive, though 
inconclusive, results [12| . In this paper we focus on the exact solution for the primordial scalar perturbation spectrum 
generated during inflation found in [l[, in which the effective mass of the inflaton, V"{(jji), experiences a sudden small 
change. This model assumes a much smoother potential than those which had previously been considered, which dealt 
either with a sudden step in the potential [^, M, \l^, or in its first derivative [^, |T0|- As a result, the corresponding 
feature in the perturbation spectrum is the weakest (in particular, superimposed oscillations are very small). Still it 
is observable, and can even produce the dominant contribution to the running of the slope of the power spectrum over 
the observable range of scales. We shall show that this new universal (i.e. not dependent on details of the change) 
feature in the infiationary perturbation spectrum is in excellent agreement with recent WMAP5 results. 



[1] Wilkinson Microwave Anisotropy Probe (WMAP) 



A microscopic model producing: such feature in the inflaton potential may be chosen similar to that used in the hybrid 
inflationary scenario [l3[ , but with different numerical values of parameters in order that a fast (quasi-equilibrium) 
second order phase transition experienced by another heavy scalar field should occur during inflation, and not at its 
end. Also this second field should be sufficiently weakly coupled to the inflaton to produce a small change, |A7tt,| ^ H, 
in the inflaton mass where H is the Hubble parameter, H = a/a, at the moment of the transition. Thus, such a 
transition may be called a 'mini-waterfall' in analogy with the 'waterfall' transition used to end inflation in the hybrid 
inflationary scenario. The sign of this coupling may be arbitrary, in principle, but we shall show that the (more 
natural) positive coupling is favoured by the WMAP5 data since it corresponds to the negative running of the power 
spectrum. 

2. INFLATION WITH MINI- WATERFALL IS TESTED USING CMB DATA 

Let us first remind the reader of the microscopic two-field model used in 1] to produce a step in the second derivative 
of an effective inflaton potential. This is just the well-known model with a Higgs-like potential, also used in the hybrid 
inflationary scenario: 

V{^P,<j^) = ^{M^- AV^')' + Im^^' + ^0V' , (2.1) 

(as pointed above, a positive value of g^ is favoured by the data, so the notation g^ for the coupling parameter is 
justified). Close to -0 = 0, the effective mass of the field ip is 

-^ - 1^ - ,^0^ - M^ . (2.2) 

At (pc — M/g the curvature of V{ip, cj)) along the ip direction vanishes so that to?, > for (p > (j)^ while to,?, < for 
(j) < <f>c- This implies that for large values of the infiaton (j) the auxiliary field tp rolls towards %p — Q. However, once 
the value of (p falls below (j)c the ^ = 0, configuration is destabilized resulting in a rapid cascade (mini- waterfall) which 
takes ip from -0 = to its zero-temperature equilibrium value "0^ = (M^ — g'^(p'^)/X where its mass becomes positive 
and large once more. We assume that the field ip is sufficiently heavy, so it quickly relaxes to its equilibrium value 
both before and after the mini- waterfall during a characteristic time 5t <C H~^ . For this to occur, the parameters of 
the model (|2.ip should satisfy a number of inequalities which were presented and investigated in [Ij. In particular, g^ 
should be much less than A. 

If so, quantum fiuctuations of the field ip may be neglected. Then ip may be excluded from the action of the model 
(|2.ip . and its non-zero equilibrium value ip{(t>) after the transition leads to the change in the effective infiaton potential: 
V{4i,%p) -^ V{(p,'ip{4>)) = V{4') (we omit tilde in what follows). Similar to free energy during a second order phase 
transition in thermodynamic equilibrium, the potential V and its first derivative with respect to (p are continuous at 
(p = <pc, while its second derivative acquires a step at this point. This is, in particular, how it is done in the theory 
of hybrid inflation [l3| (see also [11[ for a more general case). However, in contrast to hybrid inflation, the change in 
the effective mass square is small as compared to H^, so inflation continues without a break in our scenario. 

In this case, as shown in [l|, the values of the primordial spectral indices ni and n2 of the power spectrum of scalar 
perturbations induced by quantum fluctuations occurring before and after the mini-waterfall are given by 

1 /gTOp\2 K (1 — 2k) 

"' " ^ " 2^ \~W) {1 + Kf ' 

_ 1 /5TOp\2 4 + 3k-|-2k2 
"'"^ " ~2^\~W) (1 + k)2 ' 
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The parameter combinations k = 2Xm^ /g^M^ and gmp/M which appear in (|2.3p also enter into the expression for 
the number of inflationary e-folds which take place during the post-mini-waterfall period 
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In fig. [U the pair {ni, 712} is shown for different values of A/". 
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FIG. 1: Spectral indices for perturbations generated just before (ni) and immediately after (712) the phase transition are shown 
for three values of the number of inflationary e-folds occurring during the post-mini- waterfall period: TV = 40 (red, solid), 
A/" = 50 (blue, dotted) and TV = 60 (green, dashed). 

Having demonstrated that a step in the value of n^ can be accommodated within the model (|2.ip , one still needs to 
work out the precise form of the power spectrum Pi?,(fc) oc k'^'~^ . This was accomplished in [Ij where the perturbation 
equation for gauge invariant quantities was solved exactly. We summarize some of these results below. 

In the absence of a phase transition, the power spectrum of cosmological perturbations generated during inflation 
has the form Vno{k) oc fc"^"^. As shown in [l|], the effect of the phase transition is to alter V-R.o{k) by a factor ja — /?p, 
so that the final spectrum becomes 



Pfl,(fc)ap7^o(fc)x|«-/3p 



(2.5) 



where n2 is the value of the spectral index at fc :s> /cq, fco being the location of the feature. The Bogoliubov coefficients 
a, j3 can be expressed in closed analytical form as follows [l| 



a ~ (3 ~ 
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where A = ^2 — Mi = ("^1 — n2)/2, ^1,2 = | — j§7 + 3eo, and =F denotes the value of V" = ^^ at f = to T 0, to being 
the time of the phase transition. The resulting form of the spectral index in the vicinity of the feature has a step-like 
discontinuity modulated by small oscillations, as shown in figured 



RESULTS 



Fig. [3] compares our model (local running) with a Power Law model (PL). Results shown are for the best fit values 
of the cosmological parameters determined assuming flat ACDM as the background metric. The package CosmoMC 
[14| was used to confront our model with the five-year observational data from WMAP. We find that our model gives 
a better fit with Axlft ~ —3.052 as compared to the best PL model with a constant rij. This improvement in the 
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FIG. 2: A step in the second derivative of the inflaton potential leads to a step in the spectral index as shown in this figure 
which plots the primordial spectral index Ua as a function of x — k/ko- The step in ns at a; ~ 1 is followed by oscillations with 
decreasing amplitude. The parameters shown here correspond to ni = 0.97, n^ = 0.95, which agree well with WMAP5 data. 
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FIG. 3: Comparison of our model (blue, dashed) with a pure power- law model (red, solid), for the best fit values of parameters. 
The WMAP5 binned data with related error bars are also plotted for comparison. 



value of the likelihood has been obtained by introducing two more free parameters in our model as compared to the 
PL model, therefore it is reasonable. In fig. [H the solid lines are marginalized probabilities of ni (the spectral index 
before the phase transition) and the dotted line gives the marginalized probabilities of n2 (the spectral index after 
the phase transition). Observational data for multipoles with I > 2 have been used. A contribution from primordial 



tensors, satisfying the one-inflaton consistency relation, rit — — r/8, has been taken into account, too. The pivot point 
is set at 0.05Mpc~^. 

From fig. m it is clear that a value of ni slightly larger than that without a step is favoured by the data. The 
marginalised probability for n2 peaks at 0.947, while the marginalized probability for ni peaks at 0.97. Looking at 
the marginalized probability itself, we can see why a local running is suggested by our model: the peak of ni lies 
at a value significantly larger than n2. It is also interesting that the {^1,712} pair, favoured by CMB data, can be 
accommodated by the model (|2.1[) for very reasonable values of A/" lying in the interval 40 - 60, as shown in fig. [TJ In 
particular, to get just the peak values for ni and n2 given above, one has to choose gmp/M « 0.42 and k w 3.8 that 
leads to A/" ~ 40. The distributions for (rii — 712) (red solid line) and also for (1 — (ni + n2)/2) (green dashed line) 




FIG. 4: Marginalized posterior distributions for the spectral indices ni and n2. 



marginalized over all the other parameters are given in the left panel of fig. El It is quite clear from the plot that 
the data prefers (ni — 722) be smaller than (1 — (rii + n2)/2). Interestingly, the case (rii — 722) < (1 — (rii + n2)/2) 
corresponds to k > 2 in our model. That is, the potential energy of the massive infiaton m?4>'^ /2 dominates the 
vacuum energy density V{Q, 0) = M'*/4A at the moment of the phase transition. In the right panel of fig. [5] one can 
see the distribution of [{rii — 712) — (1 — (rii + "2)/2)] marginalized over all the other parameters. It is clear that the 
peak is at a value less than zero confirming the above result. 

Note that in our previous theoretical paper |l|, we were more concerned with the case k, < 1 since it leads to a 
larger value of tt-i — 71.2, in particular, the spectrum even becomes blue-tilted for fc — > if k < 1/2. But observational 
data appear to have made their verdict for the opposite case k > 1 when the potential V{il), (f>) is dominated by the 
m?4>^ /2 term during the phase transition. 

However, all formulae in [1[ and in Sec. 2 are valid for any value of k provided it is not very large. In more detail, the 
assumption of a fast phase transition is valid if m\ ^ H^ during the transition, which requires 1 + k ^ Xmrnj^M'^ 
(Eq. (4.9) of the paper [l|). For the typical values A ^ 0.1, m ~ lO^^mp, M -- 10^^77ip (see the table 1 in [J), this 
leads to K ^ 100 and g^ <§; 10~^, while k = 1 corresponds to g^ = 2 x 10^^. It can be verified that if this condition 
is satisfied, then the contribution of tp particles created during this transition to the total energy density is small 
compared to the change in the equilibrium value of V^iplt)^ 4'i't)) due to the step in its second derivative for the time 
period 5t ^ 1/ H after the transition during which the correction (|2.5ll2.8p to the quasi-flat power spectrum ViiQ{k) is 
generated, irrespective of whether these created ip particles have decayed into other particles over this time period or 
not. 

Another inequality bounding k from above follows from the requirement that the change in the spectral index 
ni — 7i2 should be much more than next-order slow-roll corrections to rii and ri2 (not taken into account in our 
calculations), which are ^ A/"^^. This leads to k <C A/" ~ 50. So, both upper limits on k needed for the validity of 
Eq. (|2.3p are effectively the same. Actually, they exclude the region only where the spectral feature itself becomes 
too small to be observable. 

Fig. [HI gives the probability distribution for fcg (location of the feature) which is larger for smaller fco values. It is 
clear that the feature, if exists at all, should preferentially lie on large scales: the marginalized upper limit for fco is 
0.00355 Mpc"^ at 95% CL. Therefore, our model suggests that the main evidence for large running in the WMAP 
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FIG. 5: Marginalized posterior distributions for (ni — 77,2) and also for (1 — (ni + n2)/2). The right panel gives the marginalized 
distribution of [(ni — 712) — (1 — (^i + n2)/2)]. 



dataset comes from sufficiently low multipoles with / < 40.^ 




ko(Mpc-') 

FIG. 6: Marginalized posterior distributions for ko- 

Marginalized probabilities obtained for other cosmological parameters using the primordial spectrum given by our 
model are shown in fig. [7] It follows there are no significant changes in derived values of the cosmological parameters 
in comparison with the results, obtained by the WMAP team, assuming a power-law model of the primordial spectrum 
i- 

4. CONCLUSIONS AND DISCUSSION 

In this paper we confront the WMAP 5-year data with an exact solution [1] for the primordial power spectrum of 
perturbations generated when the inflaton effective potential V{(j)) has a sudden small change in its second derivative 
with respect to 0, i.e. in mi. The spectrum possesses a local universal feature having the form of a step in the 
primordial spectral index rig, modulated by comparatively weak oscillations (by universal is meant that the form of 
the spectral feature does not depend upon the structure of the discontinuity in V{(j))). It results in a large (but local) 



[2] Here we use the correspondence between k and l^fj in the same form as the WMAP team: kRh = l^fj with R^ ~ f4000 Mpc. 




0.022 0.024 





13.6 14 
Age/GYr 



60 65 70 75 80 
H. 




FIG. 7: Marginalized probabilities of cosmological parameters in our model. 



running of n^. 

The simplest microscopic realization of such behaviour of V{(l)) is provided by a model having an auxiliary heavy 
scalar field which experiences a rapid second order phase transition (a mini- waterfall) during inflation in the observably 
accessible range of scales. Coupling this field to the inflaton leads to the desired type of local discontinuity in V{(j)). 
The model is similar to that used in the hybrid inflationary scenario, but in contrast to the latter, its parameters are 
chosen in such a way that: (i) the transition occurs during inflation and not at its end, (ii) the change in the inflaton 
mass is small compared to the Hubble parameter at this moment. That is why we call it a 'mini-waterfall', in contrast 
to a 'waterfall' which provides an end to inflation in the hybrid scenario. Should such a mini-waterfall be detected 
through the corresponding feature in the primordial power spectrum of scalar perturbations, it would provide direct 
experimental evidence for the naturalness of a similar (though larger) waterfall in the hybrid scenario (the latter not 
being directly observable since it lies at a very small comoving scale) . 

We find that the best xtff ^'^^ ^^^^ model shows an improvement by 3.052 over the best fit obtained assuming a 
featureless power law for the primordial spectrum. (This improvement comes at the cost of introducing two additional 
parameters.) It is shown that such a feature in the primordial spectrum, if exists at all, should lie on large scales 
ko < 0.003 Mpc~^. Anyway, this feature is not excluded by the present observational data. Better data expected 
from future CMB experiments will help to settle the question about its existence. 

An interesting problem not considered in this paper is the amount of non-Gaussianity in the statistics of primordial 
perturbations. Here, strictly speaking, one needs to distinguish two different, though related cases: the single-field 
inflationary model with an effective inflaton potential having the studied type of local non-analytic behaviour, and the 
two-field model (j2.ip . These models produce the same results for the power spectrum in the leading approximation, 
but may become non-equivalent at the level of deviations from Gaussian behaviour. Since the second field in the 
latter model is in the fast-rolling regime, one may expect larger amount of non-Gaussianity for it. We hope to return 
to this question elsewhere. 
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